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We present the first theoretical investigation of carrier-wave Rabi flopping in real atoms by em¬ 
ploying numerical simulations of high-order harmonic generation (HHG) in alkali species. Given 
the short HHG cutoff, related to the low saturation intensity, we concentrate on the features of the 
third harmonic of sodium (Na) and potassium (K) atoms. For pulse areas of 27r and Na atoms, a 
characteristic unique peak appears, which, after analyzing the ground state population, we correlate 
with the conventional Rabi flopping. On the other hand, for larger pulse areas, carrier-wave Rabi 
flopping occurs, and is associated with a more complex structure in the third harmonic. These new 
characteristics observed in K atoms indicate the breakdown of the area theorem, as was already 
demonstrated under similar circumstances in narrow band gap semiconductors. 


It is well known that semiconductors, when modeled 
as a two-level system, develop a periodic oscillation of 
the population inversion when interacting with constant 
light, a phenomenon predicted by I. I. Rabi in the 30s, 
called Rabi flopping [l| . Rabi flopping has also been ob¬ 
served when using ultrafast optical pulses e.g. 0S- For 
these pulses, peculiar behavior emerges when the driven 
light intensity is so high that the period of one Rabi os¬ 
cillation is comparable with that of one cycle of light. 
In this case, the area theorem has been shown to break 
down Q, and a new phenomena, known as carrier-wave 
Rabi flopping (CWRF), emerges. These features can be 
schematically observed in the so-called Bloch sphere (for 
the definition and more details see e.g. 0), presented in 
Fig. 1. 

In particular, Fig. 1(a) depicts the conventional Rabi 
flopping on a Bloch sphere. For this case the Rabi pe¬ 
riod is much larger than the driven light period and the 
Bloch vector, formed by the real (it) and imaginary ( v ) 
parts of the optical polarization and the population in¬ 
version (w) of a two-level system, spirals up starting from 
the south pole (corresponding to all the electrons in the 
ground state), reaches the north pole and returns to its 
initial position for the case of squared-shaped pulses with 
an envelope area of 0 = 2tt. Here optical oscillations are 
mapped to an orbit of the Bloch vector parallel to the uv 
or equatorial plane. Additionally, oscillations of the pop¬ 
ulation inversion are given by the motion in the uw plane. 
The corresponding spectrum of the optical polarization 
would exhibit then two peaks centered around the two- 
level transition frequency. On the other hand, Fig. 1(b) 
presents results for a much shorter pulse, where the Rabi 


period is equal to the driven light period. Even when the 
envelope area for this case is 0 = 4n, it is clear that the 
Bloch vector does not return to the south pole, as may 
be expected. To the contrary, a more chaotic behavior 
is observed in the motion of the Bloch vector, resulting 
in a more complex shape in the spectrum of the optical 
polarization. Furthermore, the well-known area theorem 
of the nonlinear optics fails when this parameter regime 
is reached. Note that multipeak splitting of the reso¬ 
nance fluorescence spectrum by short pulses in the stan¬ 
dard Rabi flopping regime was predicted in Refs. 00, 
although this effect is due to a complex temporal inter¬ 
ference effect, rather than CWRF. 



FIG. 1. Sketch of the Bloch Sphere showing the different 
regimes, (a) schematic showing the travel of the Bloch vector 
for conventional Rabi flopping for a pulse with an envelope 
pulse are 0 = 2-7T. (b) same for carrier-wave Rabi flopping for 
a pulse with an envelope pulse area of 0 = 47r. The bottom 
panels show the evolution of the population inversion w (see 
the text for more details). 
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Experiments on narrow band gap semiconductors have 
shown a clear signature of CWRF, which manifested itself 
as a split in the third harmonic of the emitted light into 
the forward direction [9:]. Recently, Rabi flopping and the 
consequent coherent pulse reshaping has been experimen¬ 
tally observed in a quantum cascade laser 0 , suggesting 
a new promising approach to short pulse generation. One 
of the advantages of atoms (relative to semiconductors) 
is the possibility to employ longer laser pulses and, as a 
consequence, to explore a broader range of laser parame¬ 
ters, as well as provide an alternative to carrier envelope 
phase (CEP) characterization. In addition, it has been 
shown that in semiconductors the Coulomb interaction of 
carriers in the bands gives rise to an enhancement of the 
external laser field and consequently the envelope pulse 
area by as much as a factor of two, considerably cornpli- 
cating the observed interpretation of CWRF phenomena 

iiS 

When the atom was simplistically modeled as a two- 
level system, conventional Rabi floppin g b ehavior and 
CWRF features were observed (see e.g. ]3, H3| and ref¬ 
erences therein). However, it is well known that the 
two-level approximation breaks down when strong elec¬ 
tric fields are applied, in particular in the CWRF regime. 
An important question emerges as to what extend CWRF 
could potentially be observed in real atoms. In this Let¬ 
ter we demonstrate for the first time how the CRWF sig¬ 
natures show up in the high-order harmonic generation 
(HHG) spectra of real atoms. In particular, using a ro¬ 
bust theoretical approach that accurately models both the 
ground and excited states of K atoms combined with real¬ 
istic laser parameters, we observe clearly distinct features 
in the third harmonic and correlate it with the behavior 
of the ground state population. The latter approach is 
closely related to the description of semiconductors and 
atoms modeled as two-level systems. In order to further 
support our conclusions of CWRF-likc behavior in K, we 
also compute HHG of Na atoms, in which the transi¬ 
tion energy between the ground and first excited state 
is not resonant with the driven light, and, as a conse¬ 
quence, a conventional HHG spectra, i.e. single peaks 
at odd harmonics of the driven frequency, is observed. 
Our predictions can be tested experimentally using cur¬ 
rently available ultrashort laser pulses of a Ti:Sa laser 
with wavelengths centered in the range of 750 — 800 nm. 

We start by describing our theoretical approach, 
putting special emphasis on the choice of the atomic po¬ 
tentials, such that our results for both ground and excited 
states are in excellent agreement with experimental mea¬ 
surements (see Table I). To find CWRF signatures, we 
focus on the HHG spectra. Since the HHG spectra is 
proportional to the electron dipole moment, we can es¬ 
tablish a one to one correspondence between the media, 
modeled as a collection of oscillators, and single atoms il¬ 
luminated by a strong laser field (any macroscopic effect, 
such as phase matching, could be safely neglected, con¬ 


sidering the low-order harmonic cutoff developed in alkali 
atoms, closely related to their low saturation intensity). 
To create the conditions for CWRF, we used an atomic 
system in which the period of a Rabi oscillation 0 (cor¬ 
responding to the transition between the ground and the 
first excited states 00) is similar to one period of 
the laser light. For the usual (Ti:Sa) laser sources such 
system was provided by K atoms, with a transition en¬ 
ergy between the ground and the first excited state of 
1.61 eV (hence close to the laser source photon energy of 
1.55 eV). 

Alkali atoms, due to their atomic structure (gas no¬ 
ble structure plus only one external electron), are well- 
suited to be described by the single active electron ap¬ 
proximation (SAE). We therefore focus on such atoms 
to avoid the possible role of electron-electron correla¬ 
tions, which have been found to have an important, 
yet still poorly understood role in HHG spectra 0. 
Based on SAE approximation, we use the atomic poten¬ 
tial reported in |l7[ to describe K and Na atoms. Using 
a Hartree-Fock-based method we set the two parame¬ 
ters V c and V e in the following generic potential form, 
V K ,Na(r ) = (rH ^ o)ii - -pry^), where V c accounts for the 
effect of the atomic core (nucleus plus all complete shells), 
V e represents the external potential and ro = y L . Using 
this method we find the ground, 3s, and the first ex¬ 
cited state, 3 p, of K, as well as the 4s and 4 p for Na, 
with precision of A E « ±0.0084eV (for details see Table 
I). In addition, we also compute numerically the element 
transition dipole d ns ^ np = (ip n s\z\ipnp) for both atoms 
(n = 3 for Na and n = 4 for K) and compare them with 
the experimental values reported in Ref. 0] (see Table 
I). As can be deduced, our theoretical dipole values show 
excellent agreement with experimental measurements, in¬ 
dicating that the potential given above accurately repre¬ 
sents K and Na atoms. For finding optimal laser pa¬ 
rameters to experimentally observe CWRF in atoms, an 
accurate dipole moment is essential since Rabi flopping 
frequency is linearly proportional to the dipole moment 
as well as the electric field strength 0,0. In addition, 
HHG spectra, which is the focus of our investigation, 
is believed to be particularly sensitive to the details of 
electron dynamics inside atoms and molecules, making it 
crucial to not only accurately describe the ground state 
(as many atomic potentials in the literature do), but, in 
the case of resonant transitions due to Rabi flopping, also 
the excited state 00 - 

To compute HHG spectra, we numerically solve the 
three dimensional Time Dependent Schrodinger Equa¬ 
tion (3D-TDSE) in the length gauge and using the atomic 
potential, Vx^Nair), given above for K and Na atoms, re¬ 
spectively. The harmonic yield from a single atom is then 
proportional to the Fourier transform of the dipole accel¬ 
eration of its active electron and can then be obtained 
from the electronic wave function after time propaga¬ 
tion. Our code is based on an expansion of spherical 
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FIG. 2. 3D-TDSE harmonic spectra in K for the corresponding laser intensities I = 3.158 x 10 11 W/cm 2 (panel a), I = 
5.6144 x 10 11 W/cm 2 (panel b) and I = 1.108 x 10 12 W/cm 2 (panel c). Panels (d), (e) and (f) represent the HHG in Na for 
the same laser parameters. The insets of panels (a), (b) and (c) show a zoom of the third harmonic ui/uoo = 3 (see the text for 
details). 


Sodium 

Experimental (NIST) Numerical 
(present) 

3s 

5.139 eV 

5.135 eV 

3 P 

3.036 eV 

3.038 eV 

Transition dipole 

2.49 au 

2.40 au 

Potassium 

Experimental (NIST) Numerical 
(present) 

4 s 

4.340 eV 

4.347 eV 

4p 

2.730 eV 

2.725 eV 

Transition dipole 

2.92 au 

2.79 au 


TABLE I. Experimental and theoretical values of the energy 
gap between the ground and first excited states, joint with 
the transition dipole both for Na and K. 


harmonics, Y) m , and takes advantage of the cylindrical 
symmetry of the problem (hence only the m = 0 terms 
need to be considered). The time propagation is based 
on a Crank-Nicolson method implemented on a splitting 
of the time-evolution operator that preserves the norm 
of the electronic wave function. The coupling between 
the atom and the laser pulse in the length gauge, lin¬ 
early polarized along the z axis, is written as Vi(z,t) = 
E(t) z, where E(t) is the laser electric field defined by 
E(t) = E 0 sin 2 sin(u;ot + </>)■ E 0 is the laser electric 


field peak amplitude (£/> = \JI /To with Iq = 3.5 x 10 16 
W/cm 2 ), to o = 0.0596 a.u. (A = 765.1 nm), N the total 
number of cycles in the pulse and 4> the CEP. Further¬ 
more, T defines the laser period T = 2tt/u}q ss 2.5 fs. 
In the simulations presented here we consider the case 
N = 20, corresponding to an intensity envelope of full 
width half maximum (FWHM) of 0.36NT , (7.2 optical 
cycles^ 18 fs FWHM), and cf> = 0 (see below for details). 

As discussed in the introduction, we set input laser 
frequency the same value (in atomic units) as the energy 
corresponding to the transition 4s —> 4p of K in order 
to observe a CWRF-like behavior. In order to compare 
with a conventional situation (meaning the usual condi¬ 
tions for HHG), we use the same input laser parameters 
with Na atoms, for which the transition energy from the 
ground to the first excited state, 3s —>• 3 p, corresponds 
to 2.10 eV, and is therefore non-resonant with the laser 
frequency. 

In Fig. 2, we show the harmonic spectra computed 
from the 3D-TDSE for both K (Figs. 2(a), 2(b) and 2(c)) 
and Na (Figs. 2(d), 2(e) and 2(f)) atoms. We have cho¬ 
sen the laser parameters to cover three different regimes, 
namely, for panels 2(a) and 2(d) the envelope pulse area, 
®K,Na, is close to 27 t. The envelope pulse area is de¬ 
fined as &K,Na ~ dK,Na,EoAt , where dfc,Na is the dipole 
transition matrix for K or Na (see Table I) and At the 
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FIG. 3. Time evolution of the ground state population (red thick line) along the laser pulse (blue thin line) corresponding to 
the cases plotted in Fig. 2. 


FWHM pulse duration (At = 18 fs i.e., At « 750 au). In 
particular, for a laser intensity I = 3.158 x 10 11 W/cm 2 
(Eq = 0.003 a.u.), 0 k ~ 27t and 0 jv q ~ 5.4 (for com¬ 
parison see the values used in semiconductor GaAs [§j]). 
Panels 2(c) and 2(f) correspond to values of Qk and 
0AT a , respectively, close to 4-7T (11.7 for K and 10.1 for 
Na), obtained by keeping constant the pulse duration 
and now using a laser intensity I = 1.108 x 10 12 W/cm 2 
(Eq = 0.0056 a.u.). Panels 2(b) and 2(e) were chosen to 
have an intermediate value of intensity I = 5.6144 x 10 11 
W/cm 2 (E 0 = 0.004 a.u.), corresponding to pulse enve¬ 
lope areas of 8.4 and 7.2 for K and Na, respectively. 

From the HHG spectra of K atoms, we observe a dras¬ 
tic change around the third harmonic, u/w o = 3, as pulse 
envelope area increases. Our results are directly analo¬ 
gous to the manifestation of CWRF behavior observed 
in semiconductors in [3]. Note the marked contrast to 
the HHG spectra of Na, shown in panels 2(d), 2(e), and 
2(f), where a characteristic peak is present in the third 
harmonic regardless of the envelope pulse area. The on¬ 
set of this more complex behavior in K atoms for suffi¬ 
ciently large envelope pulse areas is in agreement with 
conclusions in Q, where the onset of CWRF behavior 
and the consequent break-down of the area theorem was 
predicted for a two-level resonant systems when the pulse 


envelope area significantly exceeds 27 t. 

To get further insight into the physical mechanism be¬ 
hind the complex structure of the HHG spectra in K 
atoms, in Fig. 3 we present the time dynamics of the 
ground state population for all the cases depicted in 
Fig. 2. For a two-level system, used as a prototypical 
model for a semiconductor or a simplistic picture for a 
real atom, the electron dynamics due to interaction with 
laser light can be represented on a Bloch sphere (for de¬ 
tails see e.g 01)- In this case, the ground state popu¬ 
lations and the regular Rabi oscillations can be depicted 
as moving along the surface of the sphere (see Fig. 1(a)). 
When the CWRF regime is reached, clear signatures, 
corresponding to the break-down of the area theorem, 
should occur on the Bloch sphere as well (see Fig. 1(b)). 

Following an analogy with a two-level system, we can 
distinctly observe CWRF-like behavior in Figs. 3(b) and 
3(c) and the corresponding counterpart in the HHG spec¬ 
tra (Figs. 2(b) and 2(c)). On the contrary, a conventional 
behavior in the third harmonic (Figs. 2(a) and 3(a)), can 
be correlated with: (i) ordinary Rabi oscillations for the 
case of K (Fig. 2(a)), i.e. the ground state is completely 
depopulated, even though the laser intensity is low and 
this would be analogous to a travel of the Bloch vector 
from the south to the north pole jsj; (h) normal behaviour 
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of atoms in strong field for all the Na cases (Figs. 2(d)- 
2(f)), i.e. gradual depopulation of the ground state due 
to laser ionization (Figs. 3(d)-3(f)). 




FIG. 4. HHG for K (panel a) and Na (panel b) for different 
CEPs. The laser parameters are the same as in Figs. 2(b) and 
2(e), respectively. Solid line 0 = 0, dotted line <j> = n/2. 

In conclusion, we find the signatures of CWRF in real 
atoms, by studying the third harmonic of alkali atomic 
species. Analogous to the case of semiconductors, we can 
correlate this new feature with the complex dynamics of 
the ground state population. Our model uses accurate 
values for the atomic wavefunction of both ground and 
excited states (as is evidenced by the excellent agreement 
between the calculated states energies with experimental 
values) as well as accurate laser parameters, easily achiev¬ 
able with the current laser technology. In particular, a 
Ti:Sa laser provides laser pulses with wavelengths cen¬ 
tered in the range 750 — 800 nm, very close to the 765 
nm value corresponding to the transition energy 4s —> 4 p 
in K. As a consequence, the experimental confirmation 
of our results appears straightforward. Moreover, the 
CWRF phenomenon in atoms could emerge as a robust 
alternative for CEP characterization for long pulses as 
can be seen in Fig. 4 where we show HHG spectra for K 
(4(a)) and Na (4(b)) and for two different values of the 
CEP. Note that in the case of Na (out of resonance) the 
third harmonic does not present appreciable differences. 
However the third harmonic of K is strongly affected, 
in spite of the fact the driving laser is rather long (20 
cycles) of total duration. It is well known that in a con¬ 
ventional situation the HHG spectra is only sensitive to 


the CEP changes when the driving laser field is a few- 
cycle pulse. Furthermore the CWRF could be used to 
control the laser-induced ionization, known as a crucial 
ingredient for harmonic propagation, via manipulation of 
the ground state population. 

We acknowledge the financial support of the MICINN 
projects (FIS2008-00784 TOQATA, FIS2008-06368-C02- 
01, and FIS2010-12834), ERC Advanced Grant QUA- 
GATUA and OSYRIS, the Alexander von Humboldt 
Foundation (M.L.), and the DFG Cluster of Excel¬ 
lence Munich Center for Advanced Photonics. This re¬ 
search has been partially supported by Fundacio Privada 
Cellex. J.A.P.-H. and L. Roso acknowledge support from 
Laserlab-Europe (Grant No. EU FP7 284464) and the 
Spanish Ministerio de Economia y Competitividad (FU- 
RIAM Project FIS2013-47741-R). We thanks Christian 
Hackenberger for helping us with the artwork. 


* marcelo.ciappina@mpq.mpg.de 

[1] I. I. Rabi, Phys. Rev. 49, 324 (1936). 

[2] S. T. Cundiff, A. Knorr, J. Feldmann, S. W. Koch, E. O. 
Gobel and H. Nickel, Phys. Rev. Lett. 73, 1178 (1994). 

[3] S. Hughes, Phys. Rev. Lett. 81, 3363 (1998). 

[4] M. Wegener, Extreme Nonlinear Optics (Springer-Verlag, 
Berlin, 2005). 

[5] K. Rzazewski and M. Florjanczyk, J. Phys. B 17, L509 
(1984). 

[6] M. Florjanczyk, K. Rzazewski and J. Zakrzewski, Phys. 
Rev. A 31, 1558 (1985). 

[7] K. Rzazewski, J. Zakrzewski, M. Lewenstein and J. W. 
Haus, Phys. Rev. A 31, 2995 (1985). 

[8] M. Lewenstein, J. Zakrzewski and K. Rzazewski, J. Opt. 
Soc. Am. B 3, 22 (1986). 

[9] O. D. Miicke, T. Tritschler, M. Wegener, U. Morgner and 
F. X. Kartner, Phys. Rev. Lett. 87, 057401 (2001). 

[10] H. Choi, V.-M. Gkortsas, L. Diehl, D. Bour, S. Corzine, 
J. Zhu, G. Holler, F. Capasso, F. X. Kartner and T. B. 
Norris, Nat. Phot. 4, 706 (2010). 

[11] C. Ciuti, C. Piermarocchi, V. Savona, P. E. Selbmann, 
P. Schwendimann and A. Quattropani, Phys. Rev. Lett. 
84, 1752 (2000). 

[12] M. Frasca, J. Opt. B 3, S15 (2001). 

[13] L. Allen and J. H. Eberly, Optical Resonance and Two- 
level Atoms (Wiley, 1975). 

[14] B. Sundaram and P. W. Milonni, Phys. Rev. A 41, 6571 
(1990). 

[15] P. Meystre, Opt. Commun. 90, 41 (1992). 

[16] A. D. Shiner, B. E. Schmidt, C. Trallero-Herrero, H. J. 
Worner, S. Patchkovskii, P.B. Corkum, J.-C. Kieffer, F. 
Legare and D.M. Villeneuve, Nat. Phys. 7, 464 (2011). 

[17] J. A. Perez-Hernandez, Fotoionizacion de atomos alcali- 
nos mediante un laser pulsado intenso (Master Thesis, 
Universidad de Salamanca, 2004). 

[18] D. A. Steck, Los Alamos National Laboratory. Los 
Alamos, NM 87545. http://steck.us/alkalidata/ 

[19] A.E. Boguslavskiy, J. Mikosch, A. Gijsbertsen, M. Span¬ 
ner, S. Patchkovskii, N. Gador, M.J..I. Vrakking and A. 
Stolow, Science 335, 1336 (2012). 

























